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Core Noise Diagnostics of Turbofan Engine Noise Using 
Correlation and Coherence Functions 

Jeffrey H. Miles 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Summary 

Cross-correlation and coherence functions are used to look for periodic acoustic com- 
ponents in turbofan engine combustor time histories, to investigate direct and indirect 
combustion noise source separation based on signal propagation time delays, and to pro- 
vide information on combustor acoustics. Using the cross-correlation function, time delays 
were identified in all cases, clearly indicating the combustor is the source of the noise. In 
addition, unfiltered and low-pass filtered at 400 Hz signals had a cross-correlation time de- 
lay near 90 ms, while the low-pass filtered at less than 400 Hz signals had a cross-correlation 
time delay longer than 90 ms. Low-pass filtering at frequencies less than 400 Hz partially 
removes the direct combustion noise signals. The remainder includes the indirect com- 
bustion noise signal, which travels more slowly because of the dependence on the entropy 
convection velocity in the combustor. Source separation of direct and indirect combustion 
noise is demonstrated by proper use of low-pass filters with the cross-correlation func- 
tion for a range of operating conditions. The results may lead to a better idea about the 
acoustics in the combustor and may help develop and validate improved reduced-order 
physics-based methods for predicting direct and indirect combustion noise. 


N omenclat ur e 

B e resolution bandwidth, Hz, B e = 1/T \ j = r s /N = 2 Hz 
c c combustor region speed of sound, m/s 

D propagation time delay or lag, s 

f c upper frequency limit, f c = 1/2 At = r s / 2, Hz (32 768 Hz) 

/l cutoff frequency of low-pass filter 

H L p[f L ](juj) low-pass filter transfer function with cutoff frequency, /l 

L c combustor region length, m 

L y number of frequencies, / C /A/ = N/2 (16 384) 

M c combustor region Mach number 

N segment length, number of data points per segment (32 768) 
n 0 number of overlapped data segments/blocks 

n s number of disjoint (independent) data segments/blocks, n s = ^total ~ 128 

Pi confidence interval percent value used to calculate coherence threshold 

r microphone radial location, m (30.48 m) 

R x autocorrelation function 

r s sample rate, samples/s (65 536) 

R xy cross-correlation function 

T^otal record length, s (« 70 s) 

t time, s 

T c temperature in combustor, °K 

Td(i) record length of segment z, AT/r s , 0.5 s 
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v s combustor region flow convection velocity, m/s 

W molecular weight of air, 28.97 kg/k mol 

x time-dependent signal 

y time-dependent signal 

Subscripts 

i running segment index 

Symbols 

A t(/l) Ipeak measured combustor convection time delay less 90 ms found using low-pass filter, s 
A f frequency step, l/T^ = r s /N , Hz (2 Hz) 

At sampling interval, 1 jr s (1/65 536), s 

At combustor region time delay, s 

7 ratio of specific heats, 1.4 

7 2 n magnitude squared coherence function of noise 

7 %y magnitude squared coherence function of signals x and y 

3? gas constant, 8314 J/k mol K 

r time, s 

r a combustor region acoustic wave travel time, s 

r s combustor region entropy travel time, s 

p xy cross-correlation coefficient 


I. Introduction 


The study of turbofan engine core noise is part of the NASA Foundational Research Program to 
develop noise diagnostic tools, noise prediction schemes, and noise reduction technologies. Combustion noise, 
both direct and indirect, is an important contributor to low-frequency core noise. Combustion noise is difficult 
to study because hydrodynamic and acoustic pressure fluctuations are intermixed at the source location. In 
addition, much of the far-field noise is not correlated with the combustor and/or turbine internal noise. 
However, the part that is correlated can be studied using signal cross-correlation functions and coherence 
functions. The goal of this study is to improve understanding of the acoustics in combustors that ultimately 
will enable development of improved reduced-order physics-based methods for predicting direct and indirect 
combustor noise. 

The report compares time domain and spectral domain methods using cross-correlation functions 
and coherence functions. The methods are applied to data from a dual-spool TECH977 turbofan engine. 
Acoustic data from the same TECH977 engine test program are discussed by Miles (2008, 2009), 1,2 Mendoza 
et al. (2008), 3 Weir and Mendoza (2008), 4 Schuster (2008), 5 Royalty and Schuster (2008), 6 Dougherty and 
Mendoza (2008), 7 and Weir (2008). 8 The research discussed is an extension of the study of the combustion 
noise of the TECH977 engine conducted by Miles (2008), 1 using a coherent output power spectral domain 
coherence function approach. 

The core noise components of the dual-spool turbofan engine were separated by Miles (2008) 1 using 
coherence functions. A source location technique was used that adjusted the time delay between the com- 
bustor pressure sensor signal and the far-field microphone signal to maximize the coherence and remove as 
much variation of the phase angle with frequency as possible. For the 130° microphone, a 90.03 ms time 
shift worked best for the frequency band from 0-200 Hz, while a 86.98 ms time shift worked best for the 
frequency band from 200-400 Hz. Hence, the 0-200 Hz band signal took more time than the 200-400 Hz 
band signal to travel the same distance. This suggests the 0-200 Hz coherent cross-spectral density band 
is partly due to indirect combustion noise attributed to entropy fluctuations, which travel at a low-flow 
velocity in the combustor until interactions with the turbine pressure gradient produce indirect combustion 
noise. The indirect noise source was studied by Pickett (1975), 9 Cumpsty and Marble (1977), 10,11 Cumpsty 
(1979), 12 and Gliebe et al. (2000). 13 The signal in the 200-400 Hz frequency band is attributed mostly 
to direct combustion noise. The method discussed herein is successful because acoustic and temperature 
fluctuations caused by hot spot convection were found by Miles (1983a) 14 to be related by a linear transfer 
function that includes a convective time delay. This experiment involved the measurements of pressure and 
temperature disturbances in a long tube connected to a combustor. This linear connection of entropy and 
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pressure fluctuations implies the direct and indirect combustion noise are dependent. 

Cross-correlation functions are discussed in standard signal processing texts such as those by Ben- 
dat and Piersol (1966), 15 Bendat and Piersol (1971), 16 Bendat and Piersol (1980), 17 Childers and Durling 
(1975), 18 Kay (1987), 19 and Stearns and David (1993). 20 Computational methods are discussed by Childers 
and Durling (1975), 18 Kay (1987), 19 and Stearns and David (1993). 20 The cross-correlation function was 
used to detect periodic signals in noise by Lee (1950). 21 White (1969) 22 used the cross-correlation function 
and filters to look for propagation time delays in a system with periodic signals. 

Cross-correlation functions were used by Karchmer and Reshotko (1976) 23 to study propagation 
time delays from a AVCO-Lycoming YF -102 turbofan engine combustor to a far-field microphone located at 
a 120° polar angle. Combustion noise below 250 Hz was attributed to an indirect combustion noise source. 
Cross-correlation functions were obtained using signals low-pass-filtered to remove frequencies above 1600 
Hz or 240 Hz. Karchmer (1977a) 24 presented additional results in his dissertation. In particular, cross- 
correlation functions were shown again using signals low-pass-filtered to remove frequencies above 240 Hz. 
Reshotko and Karchmer (1980) 25 presented similar test results for the Pratt & Whitney JT15D turbofan 
engine. Combustion noise below 250 Hz was attributed a direct combustion noise source. Cross-correlation 
functions were shown using signals low-pass-filtered to remove frequencies above 240 Hz. Mathews et al. 
(1977b) 26 studied propagation time delays from a Pratt & Whitney JT8D-109 engine combustor to the 
far field using the cross-correlation function and attributed the noise source to the combustor. The report 
attributed the noise to direct combustion noise and cited calculations that indicate indirect combustion 
noise was much lower than direct combustion noise. Using the cross-correlation function, Shivashankara 
(1978) 27 studied propagation time delays from an auxiliary power unit. Combustion noise below 400 Hz 
was attributed by Shivashankara (1978) 27 to a direct combustion noise source. Cross-correlation functions 
were shown by Shivashankara (1978) 27 using signals filtered to remove frequencies above 1000 Hz. More 
recently, Harper-Bourne et al. (2008) 28 used the cross-correlation function to study combustion noise of a 
turbofan demonstrator engine. The engine study program used internal and external sensors. The method 
revealed delay times from the internal combustor sensor to the far-field microphone, even though the signal 
correlated contained both periodic and random components. In addition, indirect combustion noise has 
been investigated numerically and with model scale experiments by Schemel et al. (2004), 29 Richter et al. 
(2005), 30 and Bake et al . 31 - 33 

The fundamental premise of this report is that filtering to remove frequencies above 400 Hz leaves 
the combined direct and indirect combustion noise unchanged in the signal. Filtering to remove frequencies 
above 200 Hz leaves the indirect combustion noise mostly unchanged. In addition, the low-pass filters are 
designed to eliminate phase angle modifications and to leave the signal delay time identified as the location 
of a peak of the correlation function of random noise at the proper time value. 

Fundamental to the interpretation of the results presented is that entropy waves or hot spots travel 
at the flow velocity, which is only a fraction of the speed of sound. According to Hill and Peterson , 34 the 
average axial velocity of the combustor reactants, v s , is of the order of 30 m/s. The speed of sound in a 
combustor, c c , is generally calculated based on the assumption that the air is a perfect gas . 34-36 

c c = yfim yw ( 1 ) 

As reported by Hill and Peterson , 34 typical turbine inlet temperatures are of the order of 1600° K. Conse- 
quently, the speed of sound, c c , is of the order of 800 m/s. Thus, the Mach number of the hot spots is about 
0.0375. This assumption is not a practical limitation to the present investigation. However, the use of the 
perfect gas assumption is subject to certain limitations as pointed out by Barrett and Suomi (1949) 37 and 
Cramer (1993), 38 such as the assumption that the incremental pressure due to the sound wave is very small 
compared to the static pressure of the medium; the frequencies are low so that 7 is independent of frequency. 

This report examines signals at the 130° microphone. Engine operating conditions of 48, 54, 60, and 
71 percent of maximum power are evaluated. Results are presented for an unfiltered condition, and after 
low-pass filtering to remove frequencies above 400 Hz, 300 Hz, 200 Hz, 150 Hz, and 100 Hz. In addition, for 
the 400 Hz and 150 Hz low-pass filtered signals, aligned coherence functions and coherent combustion power 
results are shown using a time shift of 86.98 and 90.03 ms. In the next section, the engine noise data and 
analysis are presented. 
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II. Engine Noise Data and Analysis 


A. Engine test data 

The static engine test was conducted at Honeywell’s San Tan Outdoor Acoustic test facility. A polar 
array of 32 ground-plane microphones was used. Each B & K 4134 microphone in an inverted stand with 
a 7-mm ground spacing was placed on a 100 ft radius. The microphones were positioned from 5° to 160° 
(relative to the inlet) in 5° increments. The microphone at 130° was used in this study. 

The dual-spool, turbofan engine has a direct drive, wide chord fan connected by a long shaft to 
the low-pressure turbine spool and a high-pressure compressor connected by a concentric short shaft to the 
turbine high-pressure spool. The fan diameter is 34.2 in. The combustor design is a straight-through-flow 
annular geometry with 16 fuel nozzles and 2 igniters. One of the igniters is replaced by a pressure transducer, 
which is identified herein as CIP1. The dependent source separation technique used to identify direct and 
indirect combustion noise uses data from the CIP1 transducer and a far-field microphone. 

The data acquisition system had a sampling rate of 65 536 Hz and a sampling duration of roughly 
70 s. This permitted data reduction using 254 overlapped ensemble averages at a bandwidth resolution of 2 
Hz. Further signal estimation parameters are shown in Table 1. 

Table 1 Spectral estimate parameters 


Parameter Value 

Segment length, (data points per segment), N 32 768 

Sample rate, r s , samples/s 65 536 

Segment length, Tj = N/r s , s 0.500 

Sampling interval, At = l/r s , s 1/65 536 

Bandwidth resolution, B e = A / = l/Tj = r s /N , Hz 2.0 

Upper frequency limit, f c = l/2At = r s / 2, Hz 32 768 

Number of frequencies, Ly = / C /A/ = N/2 16 384 

Propagation time delay/lag (T = 9° C, r = 30.48 m) D = 5 927/65 536, s 0.09044 
Number of independent samples, n s 128 

Overlap 0.50 

Sample length, s ^ 70 


B. Correlation functions 

The autocorrelation function describes the general dependence of the values of data at one time on 
the value at another time. The autocorrelation function is 

1 r T 

Rxx(t) = lim — / x(t)x(t + r)dt (2) 

T->oo 1 J o 

The cross-correlation function describes the general dependence of the values of one set of data at 
one time on the value of another set of data at another time. The cross-correlation function is 


1 r 

R xy(r) = Ryx(-r) = lim — / x(t)y(t + r)dt (3) 

T->oo 1 J Q 

Note, the cross-correlation function is not symmetric in time about zero. 

The normalized cross-correlation function known as the cross-correlation function coefficient (nor- 
malized cross-covariance function) used herein to plot the cross-correlation functions is defined as 


Pxy( T ) ~ Pyx ( T ) ~ 


Rxy( r ) 

y/Rxx{Q)Ryy(fy 


(4) 


by Bendat and Piersol. 16, 17 While the total cross-correlation functions are calculated, only the positive time 
delay is shown. 
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C. Filtering 

The 100, 200, 300, and 400 Hz low-pass filters used are three-section Butterworth filters. Filter 
design was done using bilinear transformation following a method described by Stearns and David. 20 These 
filter parameters are given in Tables 2a, 2c, 2d, and 2e. 

To see if the results depended on the filter design approach, a non-Butterworth filter was used as 
a 150 Hz filter with coefficients shown in Table 2b. For the first section, this filter uses the second-section 
filter values of the 300 Hz filter shown in Table 2d. For the second and third section, this filter uses the 
third-section filter values of the 300 Hz filter shown in Table 2d. The 150 Hz filter has a larger transition 
band than a properly designed Butterworth filter. 

To make certain the filters introduce no phase angle modification, the filter for each segment is 
applied to the forward time series and then again, with the time reversed. The time series reversal is then 
removed before processing the next segment or quitting. This procedure is discussed by Kormylo and Jain 39 
and by Hamming. 40 Note the effective filter function is the square of the absolute value of filter function 
designed. 


Table 2a 100 Hz lowpass filter parameters 


Parameter Values 

Mz) 

-1.994957877158 

0.995049566793 


Bi(z) 

0.000022922409 

0.000045844817 

0.000022922409 

A 2 (z) 

-1.986441602897 

0.986532901119 


B 2 (z) 

0.000022824555 

0.000045649111 

0.000022824555 

A 3 (z) 

-1.981557756270 

0.981648830027 


B 3 (z) 

0.000022768439 

0.000045536878 

0.000022768439 


Table 2b 150 Hz lowpass filter parameters 


Parameter Values 


M{z) 

-1.959329632211 

0.960140349842 


Bi(z) 

0.000202679408 

0.000405358816 

0.000202679408 

A 2 (z) 

-1.945140180302 

0.945945026722 


B 2 (z) 

0.0002012116a05 

0.000402423210 

0.000201211605 

A 3 (z) 

-1.945140180302 

0.945945026722 


Bs(z) 

0.000201211605 

0.000402423210 

0.000201211605 


Table 2c 200 Hz lowpass filter parameters 


Parameter Values 


Ai(z) 

-1.989758186880 

0.990124031531 


Bi(z) 

0.000091461163 

0.000182922325 

0.000091461163 

A 2 (z) 

-1.972884422898 

0.973247165073 


B 2 (z) 

0.000090685544 

0.000181371087 

0.000090685544 

A 3 {z) 

-1.963272029823 

0.963633004625 


Bs(z) 

0.000090243701 

0.000180487401 

0.000090243701 


D. White noise results 

To demonstrate the filtering process, a study was conducted using white noise generated by use of 
a “Mersenne Twister” random number generator. 41 Figure 1 shows the effect of each filter on the signal 
amplitude. 
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Table 2d 300 Hz lowpass filter parameters 


Parameter Values 

M*) 

-1.984402564179 

0.985223656312 


Bi(z) 

0.0002052730330 

0.000410546067 

0.000205273033 

A 2 (z) 

-1.959329632211 

0.960140349842 


B 2 (z) 

0.000202679408 

0.000405358816 

0.000202679408 

A 3 (z) 

-1.945140180302 

0.945945026722 


Bs(z) 

0.000201211605 

0.000402423210 

0.000201211605 


Table 2e 400 Hz lowpass filter parameters 
Parameter Values 


M{z) 

-1.978892639460 

0.980348696212 


Bi(z) 

0.000364014188 

0.000728028376 

0.000364014188 

A 2 {z) 

-1.945778355457 

0.947210046928 


B 2 (z) 

0.000357922868 

0.000715845735 

0.000357922868 

A 3 {z) 

-1.927159619380 

0.928577611302 


Bs(z) 

0.000354497980 

0.000708995961 

0.000354497980 


Figure 2 shows cross-correlation functions generated using random-number-generated white noise 
and random-number-generated white noise displaced by D T = 86.98 ms. Each plot has the cross-correlation 
function of the unfiltered signals and the signals filtered by the 100 Hz low-pass filter. Figures 2 (a) through 
(d) have the results of the 400, 300, 200, and 150 Hz low-pass filters, respectively. The location of the 
peak value is constant, independent of the filter design. The width of the peak increases as the filter design 
frequency is reduced. 

E. Coherence threshold 

A minimum threshold of the magnitude squared coherence, 7 2 n was obtained using the Pi = percent 
confidence interval as discussed by Carter, 42 ’ 43 Halliday, 44 and Brillinger. 45 This threshold is given by 

7l = l-(l-^) 1/(n ‘- 1) (5) 

and represents the coherence of two independent time series, calculated using an average of n s independent 
finite length sample segments. 

Miles 46 shows that instead of relying on the confidence interval given by Eq. (5), which is based 
on a statistical theory, to obtain a threshold value for 7^ n (/), one can use a deliberately unaligned time 
history to create the threshold value. If one of the time histories is shifted by a time delay more than 
the segment /block length, Td = N/r s , then the two time histories are totally independent unless tones are 
present. This deliberate decorrelation establishes a coherence threshold and also identifies any tones in the 
signals. Shifting the signals by this time delay removes the coherence of random noise (identified frequently 
as broadband noise) but leaves the coherence of periodic functions (identified frequently as tones). Miles 46 
found out that replacing n s with n 0 in Eq. (5), where n 0 is the number of overlapping segments, gives a 
good estimate of the coherence squared noise threshold, 7 2 n (/). Thus, for a 95 percent confidence interval, 
the noise floor is 7 2 n = 0.012 (n G = 254), where the number of overlapped averages is used, rather than the 
number of independent averages. 


III. Results 

A. Cross-correlation with and without filtering 

The sequence of low-pass filters, discussed in the previous section, is applied to the combustor 
pressure transducer signal and far-field microphone signal. The signal resulting from using the 400 Hz low- 
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pass filter contains both indirect and direct combustion noise. By applying further low-pass filters in the 
sequence, (i.e., 300 Hz, 200 Hz, 150 Hz, and 100 Hz), the resulting signal will have a decreasing amount of 
direct combustion noise. The signal resulting from using the 100 Hz low-pass filter contains mainly indirect 
combustion noise. 

The cross-correlation functions between the microphone at 130° and the combustor pressure sensor 
are presented in Figs. 3 through 10. Figures 4, 6, 8, and 10 show the cross-correlation functions in the time 
region from 80 ms to 100 ms, which contain the peak value. The left side of the peak is determined by the 
direct combustion noise, which moves at acoustic speed along its full path. The right side of the peak in the 
region is controlled by the indirect combustion noise, which initially propagates with the mean flow velocity. 

The combustion noise spectrum has a limited frequency range. For this particular engine, this limit 
is about 400 Hz. Figure 2 shows that the cross-correlation of two band limited signals is smeared in the time 
domain. Consequently, one might expect to see in Figures 4, 6, 8, and 10 a cross-correlation function with 
a one smeared peak due to direct combustion noise and one smeared peak due to indirect combustion noise. 
Examination of the cross-correlation functions shows this is the case only for the 48 percent of maximum 
power operating condition (see Figs. 3 d-f and 4). 

However, as the low-pass filter cutoff frequency is decreased, more direct combustion noise, which 
controls the left side of the peak, is removed and the temporal smearing is increased. For low-pass filters 
with a cutoff frequency less than 150 Hz one sees only a single peak (Fig. 4a). In addition, the location 
of the double-peak or single-peak maximum value shifts to larger time delays as the low-pass filter cutoff 
frequency is decreased. The no-filter case and the cases filtered with the 200-400 Hz low-pass filters show 
a first peak near 90 ms and a second smaller peak near 95 ms. The peak using the 100 Hz low-pass filter 
occurs at 92 ms. 

At the three higher operating conditions, the two-peak phenomena are not observed in cross- 
correlations made using the combustor pressure sensor with the 130° microphone; only a single peak cross- 
correlation is observed. The reason that only a single peak is detected is not clear, but it is possible that the 
relative amount of indirect combustion noise in the signal is decreased, making its detection as a discrete 
peak more difficult at these three higher power settings. Another possibility is that the ripple observed at 
the 48 percent operating condition is due to a convective entropy feedback mechanism that enhances the 
indirect combustion noise signal for this case. However, in all cases, the single peak or the left peak, where 
a double peak occurs, shifts to the right as the low-pass filter cutoff frequency is decreased. 

The resulting delay time values of the peak are shown in Table 3. For the cases where no filter is 
used or the case when the 400 Hz low-pass filter was used, this time corresponds to about 90 ms. Peak values 
of the normalized cross-correlation function are shown in Table 4. 

Table 3. Time delays at peak of cross correlation of 130° microphone signal and CIP1 combustor 
pressure signal, ms 


Engine power 
level 

No 

filter 

Low-pass 
400 Hz 

Low-pass 
300 Hz 

Low-pass 
200 Hz 

Low-pass 
«150 Hz 

Low-pass 
100 Hz 

48 percent 

90.18 

90.07 

90.33 

90.76 

91.34 

92.06 

54 percent 

89.74 

90.19 

90.52 

90.90 

91.16 

92.12 

60 percent 

90.13 

90.07 

90.38 

90.71 

90.85 

91.46 

71 percent 

90.62 

90.67 

90.91 

91.08 

91.25 

92.16 


Table 4. Peak values of normalized cross-correlation function of 130° microphone signal and 
CIP1 combustor pressure signal 


Engine power 
level 

No 

filter 

Low-pass 
400 Hz 

Low-pass 
300 Hz 

Low-pass 
200 Hz 

Low-pass 
«150 Hz 

Low-pass 
100 Hz 

48 percent 

0.026 

0.237 

0.284 

0.337 

0.379 

0.429 

54 percent 

0.023 

0.188 

0.220 

0.228 

0.182 

0.182 

60 percent 

0.016 

0.177 

0.215 

0.228 

0.178 

0.169 

71 percent 

0.013 

0.132 

0.162 

0.154 

0.109 

0.098 
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B. Tones and the coherence threshold 


The method of aligned and unaligned coherence presented by Miles 46 is extended in this report 
to include cross-correlation functions. Aligned and unaligned cross-correlation and coherence functions are 
calculated. The unaligned cross-correlation and coherence functions are calculated using signals unaligned 
by 2.74 sample segment lengths, ( D T = 2.74 Tj = 1.37 s). 

Figure 11 shows the unaligned coherence function between the microphone at 130° and the combustor 
pressure sensor, CIP1, at four power settings. The presence of tones is clearly observable. The signal 
processing is done using signals low-pass filtered at 400 Hz, unaligned by 2.74 sample segment lengths 
(D t = 2.74 Td = 1.37 s) to obtain measured coherence noise floor and to observe tones. The measured noise 
floor and the statistical noise floor are in good agreement. Tones are tabulated in Table 5. The tone at 18 
Hz occurred at all power settings. 

Figure 12 shows the corresponding unaligned cross-correlation functions. The tone at 18 Hz has a 
period of T — 1/18 = 55 ms. This wave structure is shown in Fig. 12 at each power setting. This wave is 
modulated by waves with periods of 2 to 10 ms. 

Table 5. Tones identified with 


unaligned coherence (see Fig. 11). 


Engine power 
level 

Tone 

frequency, Hz 

Coherence 

magnitude 

48 percent 

18 

0.2228 


354 

0.15 

54 percent 

18 

0.17629 


86 

0.0574 

60 percent 

18 

0.2206 


96 

0.103 


380 

0.16848 

71 percent 

18 

0785 


398 

0.2223 


C. Filtering to remove frequencies greater than 400 Hz 

Miles (2008) 1 has shown that the cross-spectrum phase angle of unfiltered signals can be used to 
separate direct and indirect combustion noise. In this section, it is shown that the method works for the 
signals that are low-pass filtered at 400 Hz. Only the results at 48 percent of maximum power are shown. 

Figure 13 shows that the phase angle in the 200-400 Hz region (direct noise source region) is flat 
if a 86.98 ms time delay is removed. Figure 14 shows that the phase angle in the 0-200 Hz region (indirect 
noise source region) is flat if a 90.03 ms time delay is removed. This yields a 3.05 ms time delay between 
the direct combustion noise acoustic signal and the indirect combustion noise acoustic signal. 

Figures 13 and 14 show the ripple in the cross-correlation function (panel a) corresponds to the peak 
in the coherence function near 36 Hz (panel b). The broadband 36 Hz peak in panel b has a coherence value 
of 0.62. The broadband 36 Hz tone might be an entropy wave feedback tone. This broadband 36 Hz tone 
does not appear in the unaligned coherence. 

The 18 Hz and the 354 Hz tones identified using the unaligned coherence (Fig. 11) are also present 
in the aligned coherence (Figs. 13 b and 14 b). The 18 Hz tone has a coherence value of 0.27. The 354 Hz 
tone might be related to a compressor disk tone. 

Ignoring the 36 Hz tone, the coherence plots in Figs. 13 b and 14 b indicate that the combustion 
noise peaks between 100 and 200 Hz. The coherent output power plots shown in Figs. 13 d and 14 d also 
indicate that the combustion noise peaks between 100 and 200 Hz. Similar results were obtained at engine 
power settings of 54, 60, and 71 percent of the maximum power setting. 
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D. Filtering to remove frequencies greater than 150 Hz 

For the 150 Hz low-pass filtered signal peak values of the cross-correlation function are shown in 
Table 3 and occur near 91 ms. These peak values are significantly different than the values obtained using 
the unfiltered signals and the signals low-pass filtered at 400 Hz which are closer to 90 ms. 

Figure 15 c shows the phase angle in the 200-400 Hz region (direct noise source region) is flat if a 
86.98 time delay is removed. Figure 16c shows the phase angle in the 0-200 Hz region (indirect noise source 
region) is flat if a 90.03 ms time delay is removed. This yields a 3.05 ms time delay between the direct 
combustion noise acoustic signal and the indirect combustion noise acoustic signal. 

Figures 15 and 16 show the ripple in the cross-correlation function (panel a) corresponding to the 
peak in the aligned coherence function near 36 Hz (panel b). The 18 Hz tone identified using the unaligned 
coherence is also present in the aligned coherence. 

Ignoring the tone, the coherence plots in Figs. 15b and 16b again indicate that the combustion noise 
peaks between 100 and 200 Hz. The coherent output power plots shown in Figs. 15d and 16d also indicate 
that the combustion noise peaks between 100 and 200 Hz. Similar results were obtained at engine power 
settings of 54, 60, and 71 percent of the maximum power setting. 

IV. Discussion 


A. Dependent Source Separation 

An interpretation of the results shown in Table 3 can be made by plotting the time delays at the 
peak of the cross correlation of the filtered 130° microphone signal and the CIP1 combustor sensor signal. 
Let 


Ar (/^) Ipeak - T (f L ) Ipeak ^0 ms (6) 

where the plotting values of the time delays are calculated relative to 90 ms. This plot is shown in Fig. 
17. As the low-pass filter design cutoff frequency, /l, is reduced, the time delays tend to increase because a 
decreasing amount of direct combustion noise is present. The reference value of 90 ms was selected since the 
results of using no filter and a 400 Hz filter are near this value as is propagation time delay shown in Table 
1. These results show low-pass filtering can be used to separate this type of dependent source. 

B. Identification of combustion noise 

The use of a sequence of low-pass filters to study the cross-correlation of the noise from the CIP1 sensor and 
the 130° far-field microphone also show that filtering can increase the correlation function significantly. The 
peak values of the correlation function shown in table 4 are plotted in Fig. 18. A low-pass filter in the range 
of 200-300 Hz works well. When a tone is present, as it is for the 48 percent of maximum power setting, the 
correlation is increasing because the tone becomes increasingly dominant. 

C. Identification of periodic components 

For an engine power setting of 48 percent maximum power, a strong periodic signal is revealed 
by the cross-correlation function calculated using the far-field microphone and the combustor sensor. This 
signal appears with or without filtering as shown in Fig. 3. The coherence shown in part (b) of Figs. 13-16 
and the coherent output power shown in part (d) of Figs. 13-16 also indicate the presence of this tone. The 
cross-correlation function indicates the tone has period of about T = 0.025 s which means the frequency 
is about f = 1/T = 40 Hz as shown. This tone corresponds to the 36 Hz peak observable in the aligned 
coherence plots shown in panel b of Figs. 13-16. The signal may be due to a feedback mechanism. At a 
higher power setting, it is not observed. 

A somewhat similar cross-correlation function was shown in Fig. 13 of a study by Reshotko and 
Karchmer 25 of combustion noise from a Pratt & Whitney JT15D at 40 percent maximum power. The cross- 
correlation function had a period of about T = 0.07 s which indicates a frequency of about f = 1/T = 142 Hz. 
This corresponded to the peak of the largest coherence function shown in Fig. 14 of the paper by Reshotko 
and Karchmer. 25 
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The cross-correlation function identifies periodic components in a manner that clearly indicates their 
significance. 

D. Source separation 

The cross-correlation function of white noise and white noise displaced by 86.98 ms as shown in Fig. 
2 produces a spike at 86.98 ms. Close examination shows some structure to the function with features that 
vary positively and negatively by 10 -4 and 10 -5 . For filtered white noise, the cross-correlation function is 
no longer sharp as shown in Fig. 2. The resolution has been decreased by filtering. The presence of tones 
creates new issues as shown in Fig. 3. 

The plots shown in Figs. 3-9 show that two separate peaks cannot always be observed in the cross- 
correlation function. Consequently, source separation of direct and indirect combustion noise cannot always 
be reliably determined directly from the cross-correlation. However, a time delay was identified in all cases 
by the cross-correlation function, clearly indicating the combustor as the source of the noise. In addition, 
using different filter cutoff frequencies clearly separates the sources. The cross-correlation function peaks at 
a time delay depending on the blend of indirect and direct combustion noise. The unfiltered signals and the 
signals low-pass filtered at 400 Hz have time delays that peak at about 90 ms. The signals low-pass filtered 
at 150 Hz have time delays that peak at about 91 ms. 

These results suggest that previous discussions of the measurements of the cross-correlation be- 
tween the combustor and a far-field microphone by Karchmer and Reshotko, 23 Karchmer, 24 Reshotko and 
Karchmer, 25 and Mathews et al. 26 did not consider that the measured time delay using filtered signals 
caused by a blend of signals from direct and indirect combustion noise. Consequently, these previous mea- 
surements might be due more to indirect combustion noise than direct combustion noise, because filtered 
cross-correlaton functions are used by Karchmer and Reshotko, 23 Karchmer, 24 Reshotko and Karchmer, 25 
and Mathews et al. 26 and some of the direct combustion noise may have been removed by filtering. 

The cross-correlation function and coherence function methods provide a procedure to detect the 
presence of coherent indirect and direct combustion noise when both are present. However, the inherent 
smearing because of filtering renders the determination of the relative contribution of indirect and direct 
combustion noise to the total combustion noise difficult. 

E. Combustor acoustic measurements 


The low-frequency acoustic cross spectrum and the acoustic cross-correlation function have been 
identified for a range of operating conditions. The circumferential acoustic combustor modes in the TECH977 
engine are discussed by Royalty and Schuster. 6 They found that most of the acoustic energy below 500 Hz 
is contained in the plane wave mode (mode order = 0). The results presented correspond to the propagation 
of plane waves from the combustor to the far field. 

F. Limitations 

The presence of tones in the signals and the decrease in resolution of the cross-correlation function 
with the filtering prevented the cross-correlation function from identifying source separation of direct and 
indirect combustion noise by finding separate time delays for each. These limitations were overcome to some 
extent by using the low-pass filter approach. However, the presence of low-frequency tones in the filter band 
remains a problem. 

The cross-spectrum phase angle approach identifies the indirect combustion noise propagation time 
to be of the order of 90 ms, while the correlation function approach identifies the indirect combustion noise 
propagation time to be of the order of 91 to 92 ms. The cross-spectrum phase angle approach identifies the 
direct combustion noise propagation time to be of the order of 87 ms, while the correlation function approach 
identifies the direct combustion noise propagation speed to be less than 90 ms. 

The lowest cutoff frequency filter removes much of the random, broadband noise. However, this 
creates signals with a correlation function strongly influenced by the remaining tones. This makes it difficult 
to obtain a valid propagation time delay for the indirect combustion noise in the 0-200 Hz frequency band. 

A 200-400 Hz bandpass filter was used in hopes of identifying a better time delay for the direct 
combustion noise in this frequency band. However, the resulting correlation function did not have a single 
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identifiable peak. In place of the peak was a modulated waveform with many peaks that appeared in the 
80-100 ms region with an envelope that had a single peak in the region. Discussing the use of sophisticated 
notch and bandpass filters to analyze this data is beyond the scope of this report. 

The correlation function procedure currently supplies support for the idea that the indirect and 
direct core noise can be separated by the time delay due to entropy wave (hot spot) convection. However, 
the method for this data set does not produce the same result as that given by examination of the cross- 
spectrum phase angle. In addition, it is more complicated to apply. The best part about the correlation 
function procedure is that it gives definite evidence that combustion noise is part of the far-field signal. 

G. Delay time consistency 

Table 3 shows the delay times measured are more dependent on the filtering than on the operating 
condition. This can be understood in general terms using information from papers on gas-turbine engine 
instability such as those in a volume edited by Lieuwen and Yang, 36 with information on gas-turbine combus- 
tion as discussed for by Lefebvre, 47 and with information presented by Hill and Peterson 34 on the mechanics 
and thermodynamics of propulsion. 

It is assumed that the delay time only depends on the difference between the travel time of the 
entropy to the turbine and the acoustic travel time to the turbine. From the turbine to the far field, the 
signals are acoustic and travel the same path to the far-field microphone. In calculating the noise produced 
by a hot spot, one might have to consider which blade produces the noise and the further propagation of a 
hot spot. These complications will be ignored. 

The delay time can be estimated by 


At = t s - r 0 

L c jv s L c /c c 

= (L c /c c )(l/M c -l) (7) 

= (L c /c c )(l il^/ c )/il^/ c ^ L c fv s 

where M c = v s /c c is the combustor Mach number, and c c = WTJW is the local speed of sound in the 
combustor. 

The difference in the measured time delay for the direct combustion noise and the indirect combustion 
noise, At, has been shown to be about 3 ms, using the method based on the cross-spectrum phase angle. The 
use of the correlation function approach yielded a mixed value of At of the order of 2 to 5 ms. Reasonable 
values of L c , T c , and M c yield a value of 2 to 5 ms for At. However, this must be left as an exercise for the 
reader. 

The interesting point about Eq. (7) is that the delay time is mainly established by the combustor 
flow velocity and combustor length. Consequently, whatever the actual combustor length and whatever the 
combustor Mach number is, over a range of engine operating conditions, the time delay measured by the 
cross-correlation function should be fairly constant. 

V. Concluding Remarks 

Unfiltered and low-pass filtered cross-correlation functions between a combustor pressure sensor and 
a far-field microphone at 130° at the 54, 60, and 71 percent of maximum power settings showed a single time 
delay peak. However, at the 48 percent of maximum power setting a single or double peak was observed 
depending on the low-pass filter cutoff frequency. The 48 percent operating condition cross-correlation 
functions showed a double peak when using signals unfiltered and filtered with cutoff frequencies of 400, 
300, and 200 Hz. For the unfiltered cross-correlation and the cross-correlation filtered using a 400 Hz cutoff 
frequency, at each power setting the largest peak is at an acoustic time delay near 90 ms. For the double 
peak case at the 48 percent operating condition, the second peak is near 95 ms. This indicates the left side 
of the double peak or single peak is due to direct combustion noise, which travels at acoustic speeds over its 
total path length; the right side of the double peak or single peak is due to indirect combustion noise. 

Low-pass filtering with a 100 Hz cutoff frequency to remove a large amount of the direct combustion 
noise leaves the indirect combustion noise which has approximately a 92 ms time delay. This time delay is 


NASA/TM— 2009-21 5609 


11 



longer than the time delay observed with no filtering or low-pass filtering at 400 Hz, which produced a time 
delay near 90 ms. This shows some of the direct combustion noise, which is always traveling at the speed 
of sound, has been removed, leaving the indirect combustion noise which spends time moving convectively 
in the combustor at a much lower velocity than the speed of sound. As a consequence of removing part of 
the combustion noise that travels fastest, the delay time of the combustion noise remaining in the signal is 
greater. The cross-correlation function, when used with low-pass filtering, is very useful in separating the 
indirect and direct combustion noise source in this system. 

The method of aligned and unaligned cross-correlation and coherence validated the statistical coher- 
ence threshold by determining by measurement the coherence of two measured time series. After delaying 
one of them for a time delay greater than the segment or block length used in the periodogram averaging, 
the coherence threshold is found. The use of the unaligned coherence removes the random noise coherence, 
leaving the coherence of any periodic process and the noise floor. In addition, the unaligned cross-correlation 
function was calculated to view the remaining periodic signals in the time domain. 

Coherence analysis using the cross-spectrum phase angle method separates direct and indirect com- 
bustion noise using filtered or unfiltered signals. This is done by an examination of the effect that removal of 
different time delays has on the phase angle of the cross-spectrum calculated using the combustor pressure 
sensor signal and a far-field microphone signal. For the cross-spectrum between a combustor pressure sensor 
and a far-field microphone at 130°, this method indicated the direct signal had a time delay of 87 ms, and 
the indirect signal had a time delay of 90 ms. No explanation is available for this discrepancy between the 
two methods. However, the time difference for either method ranges from 2 to 5 ms. This is the expected 
time range based on available information for combustor lengths and convective flow velocities. 
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Frequency, Hz 


Figure 1. Low-pass filters acting on random noise. 
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Time, x, s 


(a) 100 and 400 Hz low-pass filter. 



Time, x, s 


(c) 100 and 200 Hz low-pass filter. 


Figure 2. Cross-correllation with no filter and 

5 700/65 536 = 86.9751 ms. 
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Time, x, s 


(b) 100 and 300 Hz low-pass filter. 



Time, x, s 


(d) 100 and 150 Hz low-pass filter. 


low-pass filters of white noise displaced by Dt 








(a) 100 Hz low-pass filter. 


(b) 150 Hz low-pass filter. 


(c) 200 Hz low-pass filter. 





(d) 300 Hz low-pass filter. 


(e) 400 Hz low-pass filter. 


(f) No filter. 


Figure 3. Normalized cross-correlation function for turbofan engine operating at a setting of 48 percent 
maximum power. 
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Time, x , s 


(a) 100-400 Hz low-pass filters. 



(b) No filter. 


Figure 4. Normalized cross-correlation calculated at setting of 48 percent maximum power in region of peak 
values. 
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(a) 100 Hz low- pass Filter 


(b) 150 Hz low-pass Filter 


(c) 200 Hz low-pass Filter 



(d) 300 Hz low-pass Filter 


(e) 400 Hz low-pass Filter 


(f) No filter 


Figure 5. Normalized cross-correlation function for turbofan engine operating at a setting of 54 percent 
maximum power. 
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(a) 100-400 Hz low-pass filters 
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(b) No filter 

Figure 6. Normalized cross-correlation calculated at a setting of 54 percent maximum power in region of peak 
values. 
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Time, x, s Time, x, s Time, x, s 


(a) 100 Hz low- pass Filter (b) 150 Hz low-pass Filter (c) 200 Hz low- pass Filter 



Time, x, s Time, x, s Time, x, s 


(d) 300 Hz low-pass Filter (e) 400 Hz low-pass Filter (f) No filter 

Figure 7. Normalized cross-correlation function for turbofan engine operating at a setting of 60 percent 
maximum power. 
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(a) 100-400 Hz low-pass filters 



(b) No filter 


Figure 8. Normalized cross-correlation calculated at setting of 60 percent maximum power in region of peak 
values. 
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(a) 100 Hz low-pass filter 


(b) 150 Hz low-pass filter 


(c) 200 Hz low-pass filter 





(d) 300 Hz low-pass filter 


(e) 400 Hz low-pass filter 


(f) No filter 


Figure 9. Normalized cross-correlation function for turbofan engine operating at a setting of 71 percent 
maximum power. 
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Figure 10. Normalized cross-correlation calculated at setting of 71 percent maximum power in region of peak 
values. 
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FREQUENCY, Hz FREQUENCY, Hz 


(a) 48 percent of maximum power. 


(b) 54 percent of maximum power. 



FREQUENCY, Hz 



FREQUENCY, Hz 


(c) 60 percent of maximum power. 


(d) 71 percent of maximum power. 


Figure 11. Deliberately unaligned coherence between microphone at 130° and CIP1 combustor pressure sensor 
using signals low-pass filtered at 400 Hz unaligned by 2.74 sample segment lengths, ( Dt = 2.74 = 1.37 s) to 
measure coherence noise floor and to observe tones. 
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(a) Cross-correlation function between microphone at 
130° and CIP1 combustor pressure sensor. 



(c) Aligned unwrapped cross-spectrum phase angle be- 
tween microphone at 130° and CIP1 combustor pressure 
sensor. 


n 



Frequency, Hz 


(b) Coherence function between microphone at 130° and 
CIP1 combustor pressure sensor. 



(d) Coherent output power at microphone at 130° us- 
ing CIP1 combustor pressure sensor and a 6 dB ground 
reflection correction. 


Figure 13. Turbofan engine power setting of 48 percent maximum power, microphone at 130°. Signal processing 
uses aligned signals filtered at 400 Hz ( D T = 5 700/65 536 = 86.98 ms). 
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(a) Cross-correlation function between microphone at 
130° and CIP1 combustor pressure sensor. 



(c) Aligned unwrapped cross-spectrum phase angle be- 
tween microphone at 130° and CIP1 combustor pressure 
sensor. 


Pw 



Frequency, Hz 


(b) Coherence function between microphone at 130° and 
CIP1 combustor pressure sensor. 



(d) Coherent output power at microphone at 130° us- 
ing CIP1 combustor pressure sensor and a 6 dB ground 
reflection correction. 


Figure 14. Turbofan engine power setting of 48 percent of maximum power. Microphone at 130°. Signal 
processing done using aligned signals filtered at 400 Hz ( D T = 5 900/65 536 = 90.03 ms). 
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(a) Cross-correlation function between microphone at 
130° and CIP1 combustor pressure sensor. 



(c) Aligned unwrapped cross-spectrum phase angle be- 
tween microphone at 130° and CIP1 combustor pressure 
sensor. 


n 



Frequency, Hz 


(b) Coherence function between microphone at 130° and 
CIP1 combustor pressure sensor. 



(d) Coherent output power at microphone at 130° us- 
ing CIP1 combustor pressure sensor and a 6 dB ground 
reflection correction. 


Figure 15. Turbofan engine power setting of 48 percent maximum power. Microphone at 130°. Signal 
processing done using aligned signals filtered at 150 Hz ( Dt = 5 700/65 536 = 86.98 ms). 
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(a) Cross-correlation function between microphone at 
130° and CIP1 combustor pressure sensor. 



(c) Aligned unwrapped cross-spectrum phase angle be- 
tween microphone at 130° and CIP1 combustor pressure 
sensor. 


n _ 



Frequency, Hz 


(b) Coherence function between microphone at 130° and 
CIP1 combustor pressure sensor. 



(d) Coherent output power at microphone at 130° us- 
ing CIP1 combustor pressure sensor and a 6dB ground 
reflection correction. 


Figure 16. Turbofan engine power setting of 48 percent of maximum power. Microphone at 130°. Signal 
processing done using aligned signals filtered at 150 Hz ( Dt = 5 900/65 536 = 90.03 ms). 
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Figure 17. Time delay at peak value of cross-correlation function of low-pass filtered signals as a function of 
filter design frequency less 90 ms. 
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Figure 18. Peak values of low-pass filtered cross-correlation function as a function of filter design frequency. 
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